Metabolic engineering with adaptive laboratory evolution for phenylalanine production by Corynebacterium glutamicum.

The mutants resistant to a phenylalanine analog, 4-fluorophenylalanine (4FP), were obtained for metabolic engineering of Corynebacterium glutamicum for producing aromatic amino acids synthesized through the shikimate pathway by adaptive laboratory evolution. Culture experiments of the C. glutamicum strains which carry the mutations found in the open reading frame from the 4FP-resistant mutants revealed that the mutations in the open reading frames of aroG (NCgl2098), pheA (NCgl2799) and aroP (NCgl1062) encoding 3-deoxy-d-arabino-heptulosonate-7-phosphate, prephenate dehydratase, and aromatic amino acid transporter are responsible for 4FP resistance and higher concentration of aromatic amino acids in their culture supernatants in the 4FP-resistant strains. It was expected that aroG and pheA mutations would release feedback inhibition of the enzymes involved in the shikimate pathway by phenylalanine and that aroP mutations would prevent intracellular uptake of aromatic amino acids. Therefore, we conducted metabolic engineering of the C. glutamicum wild-type strain for aromatic amino acid production and found that phenylalanine production at 6.11 ± 0.08 g L-1 was achieved by overexpressing the mutant pheA and aroG genes from the 4FP-resistant mutants and deleting aroP gene. This study demonstrates that adaptive laboratory evolution is an effective way to obtain useful mutant genes related to production of target material and to establish metabolic engineering strategies.

Isolation and whole-genome sequencing analysis of Escherichia fergusonii harboring a heat-labile enterotoxin gene from retail chicken meat in Okinawa, Japan.

This study aimed to reveal the prevalence of heat-labile enterotoxin (LT) gene-positive Escherichia fergusonii in retail chicken meat and genetically characterize these strains. E. fergusonii harboring LT gene was isolated from 6 out of 60 (10%) retail chicken samples in Okinawa, Japan. Whole-genome sequencing analysis revealed that LT gene-positive E. fergusonii from chicken meat and feces contain an IncFII plasmid harboring elt1AB, and suggested to spread clonally to retail chicken through fecal contamination. Additionally, it was found that these strains harbor multidrug-resistant genes on their plasmids. Their pathogenicity and continuous monitoring are required for confirmation.

Targeted A-to-G base editing in the organellar genomes of Arabidopsis with monomeric programmable deaminases.

Plastids and mitochondria are 2 intracellular organelles containing DNA-encoding partial but essential components for their roles, photosynthesis, and respiration. Precise base editing in both plastid and mitochondrial genomes would benefit their gene functional analysis and crop breeding. Targeted base editing in organellar genomes relies on a protein-based genome-editing system that uses the TALE-DNA recognition motif with deaminases. This is because the efficient delivery of guide RNA for clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 systems into organelles is currently impossible. Since TALE-based base editors used in organellar genomes are usually dimeric types, in this study, we used targeted A-to-G base editing in Arabidopsis (Arabidopsis thaliana) plastid and mitochondrial genomes with monomeric TALE-based deaminase for easier assembling of vectors. As a result, inheritable targeted A-to-G base editing of adenosine triphosphatase subunit 6-2 (atp6-2) in plant mitochondrial genomes and of 16S ribosomal RNA (16S rRNA) in plastid genomes of Arabidopsis was successfully induced by monomeric TALE-based adenine deaminase (AD) without off-target mutations. The monomeric TALE-based adenine deaminases also demonstrated a preference for editing the 8th T on the same strand from the recognition end. Phenotypic analysis showed that A-to-G conversion at 1139A of plastid 16S rRNA conferred substantial spectinomycin resistance in Arabidopsis, but not the other 2 potential-resistant mutations at 1131T and 1137T, predicted from the previous bacterial data. Our study demonstrated the feasibility of monomeric TALE-based ADs in plant organelles and their potential contribution to the functional analyses of plant organelles with easier assembling.

Chromosomal-level assembly of Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki genomes.

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.

HlyF, an underestimated virulence factor of uropathogenic Escherichia coli.

OBJECTIVES: Urinary tract infections (UTIs) are primarily caused by uropathogenic Escherichia coli (UPEC). This study aims to elucidate the role of the virulence factor HlyF in the epidemiology and pathophysiology of UTIs and investigate the dissemination of plasmids carrying the hlyF gene. METHODS: An epidemiological analysis was conducted on a representative collection of 225 UPEC strains isolated from community-acquired infections. Selected hlyF+ strains were fully sequenced using a combination of Illumina and Nanopore technologies. To investigate the impact of HlyF, a murine model of UTI was utilized to compare clinical signs, bacterial loads in the bladder, kidney, and spleen, onset of bacteraemia, and inflammation through cytokine quantification among wild-type hlyF+ strains, isogenic mutants, and complemented mutants. RESULTS: Our findings demonstrate that 20% of UPEC encode the HlyF protein. These hlyF+ UPEC strains exhibited enhanced virulence, frequently leading to pyelonephritis accompanied by bloodstream infections. Unlike typical UPEC strains, hlyF+ UPEC strains demonstrate a broader phylogroup distribution and possess a unique array of virulence factors and antimicrobial resistance genes, primarily carried by ColV-like plasmids. In the murine UTI model, expression of HlyF was linked to the UPECs' capacity to induce urosepsis and elicit an exacerbated inflammatory response, setting them apart from typical UPEC strains. DISCUSSION: Overall, our results strongly support the notion that HlyF serves as a significant virulence factor for UPECs, and the dissemination of ColV-like plasmids encoding HlyF warrants further investigation.

Isolation and Genomic Characterization of a Heat-Labile Enterotoxin 1-Producing Escherichia fergusonii Strain from a Human.

Escherichia fergusonii strains have been isolated from patients with diarrhea, but their virulence determinant has not been well elucidated. Here, we report the first isolation of a heat-labile enterotoxin 1 (LT1)-producing E. fergusonii strain (strain 30038) from a patient in Japan. The complete genome sequence of strain 30038 was determined and subjected to comparative genomics and phylogenetic analyses with 195 publicly available genomes of E. fergusonii. In addition to strain 30038, the elt1 gene was also identified in an E. fergusonii strain that is phylogenetically distinct and which was isolated from poultry in the United Kingdom. Fine genomic comparison revealed that these two strains share comparable elt1-bearing plasmids. However, an intriguing distinction arises in strain 30038, wherein the plasmid has integrated into the chromosome via a recombination process mediated by an insertion sequence. The production of active LT1 toxin by strain 30038 was verified through an in vitro assay using cultured cells. A large plasmid carrying 11 antimicrobial resistance genes was also identified in strain 30038. Our results indicate that extensive surveillance of elt1-positive E. fergusonii strains as diarrheagenic pathogens is needed. IMPORTANCE Escherichia fergusonii, a species closely related to Escherichia coli, is known to cause sporadic conditions in humans, including diarrhea. However, the critical virulence factors in E. fergusonii clinical isolates remain to be identified. This study shows the first isolation of an E. fergusonii strain carrying the elt1 gene, which encodes heat-labile enterotoxin 1, from a patient with diarrhea. Our analysis of public databases also revealed the presence of elt1-positive E. fergusonii strains isolated from poultry in the United Kingdom. Interestingly, while the elt1 gene in the poultry isolate was present on a large plasmid, in the human isolate it was integrated into the chromosome, which may confer stability on the elt1-carrying genetic element. Our findings highlight the need for extensive surveillance of elt1-positive E. fergusonii strains in livestock animals.

GINGER: an integrated method for high-accuracy prediction of gene structure in higher eukaryotes at the gene and exon level.

The prediction of gene structure within the genome sequence is the starting point of genome analysis, and its accuracy has a significant impact on the quality of subsequent analyses. Gene structure prediction is roughly divided into RNA-Seq-based methods, ab initio-based methods, homology-based methods, and the integration of individual prediction methods. Integrated methods are mainstream in recent genome projects because they improve prediction accuracy by combining or taking the best individual prediction findings; however, adequate prediction accuracy for eukaryotic species has not yet been achieved. Therefore, we developed an integrated tool, GINGER, that solves various issues related to gene structure prediction in higher eukaryotes. By handling artefacts in alignments of RNA and protein sequences, reconstructing gene structures via dynamic programming with appropriately weighted and scored exon/intron/intergenic regions, and applying different prediction processes and filtering criteria to multi-exon and single-exon genes, we achieved a significant improvement in accuracy compared to the existing integration methods. The feature of GINGER is its high prediction accuracy at the gene and exon levels, which is pronounced for species with more complex gene architectures. GINGER is implemented using Nextflow, which allows for the efficient and effective use of computing resources.

Characterization and development of a plastid genome base editor, ptpTALECD.

The modification of photosynthesis-related genes in plastid genomes may improve crop yields. Recently, we reported that a plastid-targeting base editor named ptpTALECD, in which a cytidine deaminase DddA functions as the catalytic domain, can homoplasmically substitute a targeted C to T in plastid genomes of Arabidopsis thaliana. However, some target Cs were not substituted. In addition, although ptpTALECD could substitute Cs on the 3' side of T and A, it was unclear whether it could also substitute Cs on the 3' side of G and C. In this study, we identified the preferential positions of the substituted Cs in ptpTALECD-targeting sequences in the Arabidopsis plastid genome. We also found that ptpTALECD could substitute Cs on the 3' side of all four bases in plastid genomes of Arabidopsis. More recently, a base editor containing an improved version of DddA (DddA11) was reported to substitute Cs more efficiently, and to substitute Cs on the 3' side of more varieties of bases in human mitochondrial genomes than a base editor containing DddA. Here, we also show that ptpTALECD_v2, in which a modified version of DddA11 functions as the catalytic domain, more frequently substituted Cs than ptpTALECD in the Arabidopsis plastid genome. We also found that ptpTALECD_v2 tended to substitute Cs at more positions than ptpTALECD. Our results reveal that ptpTALECD can cause a greater variety of codon changes and amino acid substitutions than previously thought, and that ptpTALECD and ptpTALECD_v2 are useful tools for the targeted base editing of plastid genomes.

Genomic characterization of endemic diarrheagenic Escherichia coli and Escherichia albertii from infants with diarrhea in Vietnam.

BACKGROUND: Diarrheagenic Escherichia coli (DEC) is a group of bacterial pathogens that causes life-threatening diarrhea in children in developing countries. However, there is limited information on the characteristics of DEC isolated from patients in these countries. A detailed genomic analysis of 61 DEC-like isolates from infants with diarrhea was performed to clarify and share the characteristics of DEC prevalent in Vietnam. PRINCIPAL FINDINGS: DEC was classified into 57 strains, including 33 enteroaggregative E. coli (EAEC) (54.1%), 20 enteropathogenic E. coli (EPEC) (32.8%), two enteroinvasive E. coli (EIEC) (3.3%), one enterotoxigenic E. coli (ETEC), and one ETEC/EIEC hybrid (1.6% each), and surprisingly into four Escherichia albertii strains (6.6%). Furthermore, several epidemic DEC clones showed an uncommon combination of pathotypes and serotypes, such as EAEC Og130:Hg27, EAEC OgGp9:Hg18, EAEC OgX13:H27, EPEC OgGp7:Hg16, and E. albertii EAOg1:HgUT. Genomic analysis also revealed the presence of various genes and mutations associated with antibiotic resistance in many isolates. Strains that demonstrate potential resistance to ciprofloxacin and ceftriaxone, drugs recommended for treating childhood diarrhea, accounted for 65.6% and 41%, respectively. SIGNIFICANCE: Our finding indicate that the routine use of these antibiotics has selected resistant DECs, resulting in a situation where these drugs do not provide in therapeutic effects for some patients. Bridging this gap requires continuous investigations and information sharing regarding the type and distribution of endemic DEC and E. albertii and their antibiotic resistance in different countries.

Escherichia cryptic clade I is an emerging source of human intestinal pathogens.

BACKGROUND: Within the genus Escherichia, several monophyletic clades other than the traditionally defined species have been identified. Of these, cryptic clade I (C-I) appears to represent a subspecies of E. coli, but due to the difficulty in distinguishing it from E. coli sensu stricto, the population structure and virulence potential of C-I are unclear. RESULTS: We defined a set of true C-I strains (n = 465), including a Shiga toxin 2a (Stx2a)-producing isolate from a patient with bloody diarrhoea identified by the retrospective analyses using a C-I-specific detection system. Through genomic analysis of 804 isolates from the cryptic clades, including these C-I strains, we revealed their global population structures and the marked accumulation of virulence genes and antimicrobial resistance genes in C-I. In particular, half of the C-I strains contained hallmark virulence genes of Stx-producing E. coli (STEC) and/or enterotoxigenic E. coli (ETEC). We also found the host-specific distributions of virulence genes, which suggests bovines as the potential source of human infections caused by STEC- and STEC/ETEC hybrid-type C-I strains, as is known in STEC. CONCLUSIONS: Our findings demonstrate the emergence of human intestinal pathogens in C-I lineage. To better understand the features of C-I strains and their infections, extensive surveillance and larger population studies of C-I strains are needed. The C-I-specific detection system developed in this study will be a powerful tool for screening and identifying C-I strains.

Helicobacter kumamotonensis sp. nov., isolated from human clinical specimens.

A Gram-stain-negative, spiral bacterium (PAGU 1991T) was isolated from the blood of a patient with diffuse large B-cell lymphoma. Phylogenetic analysis based on 16S rRNA gene sequences showed that the isolate was very closely related to Helicobacter equorum LMG 23362T (99.1 % similarity), originally isolated from a faecal sample from a healthy horse. PAGU 1991T was also very closely related to PAGU 1750 in our strain library (=CCUG 41437) with 99.7 % similarity. Additional phylogenetic analyses based on the 23S rRNA gene sequence and GyrA amino acid sequence further supported the close relationship between the two human isolates (PAGU 1991T and PAGU 1750) and the horse strain. However, a phylogenetic analysis based on 16S rRNA showed that the two human isolates formed a lineage that was distinct from the horse strain (less than 99.2 % similarity). In silico whole-genome comparisons based on digital DNA-DNA hybridization, average nucleotide identity based on blast and orthologous average nucleotide identity using usearch between the two human isolates and the type strain of H. equorum showed values of less than 52.40, 93.47, and 93.50 %, respectively, whereas those between the two human isolates were 75.8, 97.2, and 97.2 %, respectively. These data clearly demonstrated that the two human isolates formed a single species, distinct from H. equorum. Morphologically, the human isolates could be distinguished by the type of flagella; the human isolates showed a bipolar sheathed flagellum, whereas that of H. equorum was monopolar. Biochemically, the human isolate was characterized by growth at 42 °C under microaerobic conditions and nitrate reduction unability. We conclude that the two human isolates, obtained from geographically and temporally distinct sources, were a novel species, for which we propose the name Helicobacter kumamotonensis sp. nov., with the type strain PAGU 1991T (=GTC 16810T=CCUG 75774T).

Turnover of mammal sex chromosomes in the Sry-deficient Amami spiny rat is due to male-specific upregulation of Sox9.

Mammalian sex chromosomes are highly conserved, and sex is determined by SRY on the Y chromosome. Two exceptional rodent groups in which some species lack a Y chromosome and Sry offer insights into how novel sex genes can arise and replace Sry, leading to sex chromosome turnover. However, intensive study over three decades has failed to reveal the identity of novel sex genes in either of these lineages. We here report our discovery of a male-specific duplication of an enhancer of Sox9 in the Amami spiny rat Tokudaia osimensis, in which males and females have only a single X chromosome (XO/XO) and the Y chromosome and Sry are completely lost. We performed a comprehensive survey to detect sex-specific genomic regions in the spiny rat. Sex-related genomic differences were limited to a male-specific duplication of a 17-kb unit located 430 kb upstream of Sox9 on an autosome. Hi-C analysis using male spiny rat cells showed the duplicated region has potential chromatin interaction with Sox9. The duplicated unit harbored a 1,262-bp element homologous to mouse enhancer 14 (Enh14), a candidate Sox9 enhancer that is functionally redundant in mice. Transgenic reporter mice showed that the spiny rat Enh14 can function as an embryonic testis enhancer in mice. Embryonic gonads of XX mice in which Enh14 was replaced by the duplicated spiny rat Enh14 showed increased Sox9 expression and decreased Foxl2 expression. We propose that male-specific duplication of this Sox9 enhancer substituted for Sry function, defining a novel Y chromosome in the spiny rat.

Gene Recruitments and Dismissals in the Argonaut Genome Provide Insights into Pelagic Lifestyle Adaptation and Shell-like Eggcase Reacquisition.

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Evolutionary History of Sexual Differentiation Mechanism in Insects.

Alternative splicing underpins functional diversity in proteins and the complexity and diversity of eukaryotes. An example is the doublesex gene, the key transcriptional factor in arthropod sexual differentiation. doublesex is controlled by sex-specific splicing and promotes both male and female differentiation in holometabolan insects, whereas in hemimetabolan species, doublesex has sex-specific isoforms but is not required for female differentiation. How doublesex evolved to be essential for female development remains largely unknown. Here, we investigate ancestral states of doublesex using Thermobia domestica belonging to Zygentoma, the sister group of Pterygota, that is, winged insects. We find that, in T. domestica, doublesex expresses sex-specific isoforms but is only necessary for male differentiation of sexual morphology. This result supports the hypothesis that doublesex initially promoted male differentiation during insect evolution. However, T. domestica doublesex has a short female-specific region and upregulates the expression of vitellogenin homologs in females, suggesting that doublesex may already play some role in female morphogenesis of the common ancestor of Pterygota. Reconstruction of the ancestral sequence and prediction of protein structures show that the female-specific isoform of doublesex has an extended C-terminal disordered region in holometabolan insects but not in nonholometabolan species. We propose that doublesex acquired its function in female morphogenesis through a change in the protein motif structure rather than the emergence of the female-specific exon.

Targeted base editing in the mitochondrial genome of Arabidopsis thaliana.

Beyond their well-known role in respiration, mitochondria of land plants contain biologically essential and/or agriculturally important genes whose function and regulation are not fully understood. Until recently, it has been difficult to analyze these genes or, in the case of crops, to improve their functions, due to a lack of methods for stably modifying plant mitochondrial genomes. In rice, rapeseed, and Arabidopsis thaliana, mitochondria-targeting transcription activator-like effector nucleases (mitoTALENs) have recently been used to disrupt targeted genes in an inheritable and stable manner. However, this technique can also induce large deletions around the targeted sites, as well as cause ectopic homologous recombinations, which can change the sequences and gene order of mitochondrial genomes. Here, we used mitochondria-targeting TALEN-based cytidine deaminase to successfully substitute targeted C:G pairs with T:A pairs in the mitochondrial genomes of plantlets of A. thaliana without causing deletions or changes in genome structure. Expression vectors of the base editor genes were stably introduced into the nuclear genome by the easy-to-use floral dipping method. Some T1 plants had apparent homoplasmic substitutions that were stably inherited by seed progenies, independently of the inheritance of nuclear-introduced genes. As a demonstration of the method, we used it to restore the growth of an organelle transcript processing 87 (otp87) mutant that is defective in the editing of RNA transcripts of the mitochondrial atp1 gene and to identify bases in atp1 that affect the efficiency of RNA editing by OTP87.

Mining RNA-seq data reveals the massive regulon of GcvB small RNA and its physiological significance in maintaining amino acid homeostasis in Escherichia coli.

Bacterial small RNAs regulate the expression of multiple genes through imperfect base-pairing with target mRNAs mediated by RNA chaperone proteins such as Hfq. GcvB is the master sRNA regulator of amino acid metabolism and transport in a wide range of Gram-negative bacteria. Recently, independent RNA-seq approaches identified a plethora of transcripts interacting with GcvB in Escherichia coli. In this study, the compilation of RIL-seq, CLASH, and MAPS data sets allowed us to identify GcvB targets with high accuracy. We validated 21 new GcvB targets repressed at the posttranscriptional level, raising the number of direct targets to >50 genes in E. coli. Among its multiple seed sequences, GcvB utilizes either R1 or R3 to regulate most of these targets. Furthermore, we demonstrated that both R1 and R3 seed sequences are required to fully repress the expression of gdhA, cstA, and sucC genes. In contrast, the ilvLXGMEDA polycistronic mRNA is targeted by GcvB through at least four individual binding sites in the mRNA. Finally, we revealed that GcvB is involved in the susceptibility of peptidase-deficient E. coli strain (Δpeps) to Ala-Gln dipeptide by regulating both Dpp dipeptide importer and YdeE dipeptide exporter via R1 and R3 seed sequences, respectively.

Analysis of Sex Chromosome Evolution in the Clade Palaeognathae from Phased Genome Assembly.

Birds in the clade Palaeognathae, excluding Tinamiformes, have morphologically conserved karyotypes and less differentiated ZW sex chromosomes compared with those of other birds. In particular, the sex chromosomes of the ostrich and emu have exceptionally large recombining pseudoautosomal regions (PARs), whereas non-PARs are classified into two strata according to the date of their origins: stratum 0 and stratum 1 (S1). However, the construction and analysis of the genome sequences in these regions in the clade Palaeognathae can be challenging because assembling the S1 region is difficult owing to low sequence diversity between gametologs (Z-linked and W-linked sequences). We addressed this issue by applying the Platanus-allee assembler and successfully constructed the haplotype-resolved (phased) assembly for female emu, cassowary, and ostrich using only sequence read data derived from the Illumina platform. Comparative genomic and phylogenetic analyses based on assembled Z-linked and W-linked sequences confirmed that the S1 region of emu and cassowary formed in their common ancestor. Moreover, the interspersed repetitive sequence landscapes in the S1 regions of female emu showed an expansion of younger repetitive elements in the W-linked S1 region, suggesting an interruption in homologous recombination in the S1 region. These results provide novel insights into the trajectory of sex chromosome evolution in the clade Palaeognathae and suggest that the Illumina-based phased assembly method is an effective approach for elucidating the evolutionary process underlying the transition from homomorphic to differentiated sex chromosomes.

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features.

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.

Elucidation of the speciation history of three sister species of crown-of-thorns starfish (Acanthaster spp.) based on genomic analysis.

The crown-of-thorns starfish (COTS) is a coral predator that is widely distributed in Indo-Pacific Oceans. A previous phylogenetic study using partial mitochondrial sequences suggested that COTS had diverged into four distinct species, but a nuclear genome-based analysis to confirm this was not conducted. To address this, COTS species nuclear genome sequences were analysed here, sequencing Northern Indian Ocean (NIO) and Red Sea (RS) species genomes for the first time, followed by a comparative analysis with the Pacific Ocean (PO) species. Phylogenetic analysis and ADMIXTURE analysis revealed clear divergences between the three COTS species. Furthermore, within the PO species, the phylogenetic position of the Hawaiian sample was further away from the other Pacific-derived samples than expected based on the mitochondrial data, suggesting that it may be a PO subspecies. The pairwise sequentially Markovian coalescent model showed that the trajectories of the population size diverged by region during the Mid-Pleistocene transition when the sea-level was dramatically decreased, strongly suggesting that the three COTS species experienced allopatric speciation. Analysis of the orthologues indicated that there were remarkable genes with species-specific positive selection in the genomes of the PO and RS species, which suggested that there may be local adaptations in the COTS species.

Targeted base editing in the plastid genome of Arabidopsis thaliana.

Bacterial cytidine deaminase fused to the DNA binding domains of transcription activator-like effector nucleases was recently reported to transiently substitute a targeted C to a T in mitochondrial DNA of mammalian cultured cells1. We applied this system to targeted base editing in the Arabidopsis thaliana plastid genome. The targeted Cs were homoplasmically substituted to Ts in some plantlets of the T1 generation and the mutations were inherited by their offspring independently of their nuclear-introduced vectors.